Vortex Lattice Jamming and Metastability Studies by Small Angle Neutron Scattering 
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Recently extensive vortex lattice metastability was reported in MgB2 in connection with a sec- 
ond order rotational phase transition. However, the mechanism responsible for these well-ordered 
metastable vortex lattice phases is not well understood. Using small- angle neutron scattering, we 
studied the vortex lattice in MgB2 as it was driven from a metastable to the ground state through a 
series of small changes in the applied magnetic field. Our results show that metastability in MgB2 
is not the result of pinning defects. Rather, we propose that it is due to the jamming of counter 
rotated vortex lattice domains which prevents a rotation to the ground state orientation. Although 
granular jamming has been applied to the motion of individual vortices, this is the first experimental 
evidence that vortex lattice domains can act as granular entities. 

PACS numbers: 64.60. My, 74.25. Ha,74.70.Ad,61. 05. fg 



Recently, we discovered the existence of well-ordered 
metastable vortex lattice (VL) phases in MgB2 single 
crystals The VL of MgB2 consists of three hexag- 
onal phases separated by second order rotation transi- 
tions. By field cooling across the phase boundaries, it is 
possible to lock in a long lived, metastable phase. Such 
robust metastability raises the question: What mecha- 
nism is responsible for the longevity of the metastable 
states, preventing them from immediately rotating to the 
ground state? 

An obvious possible answer to the above question is 
that the metastable VLs are held in place by vortex pin- 
ning. However, pinning in connection with a VL reorien- 
tation transition was previously reported in the borocar- 
bide superconductor YM2B2C and led to a disordering 
of the VL |2j. In contrast, well ordered metastable VLs 
are seen in MgB2, making it unlikely that metastabil- 
ity is caused by pinning. Instead, it was suggested by 
Das et al. that the VL domains act as granular entities, 
jamming against one another and preventing them from 
rotating to the ground state [1]. This is distinct from 
the jamming of individual vortices observed in connec- 
tion with artificial pinning potentials Rather, it is 
analogous to jamming observed in granular systems Q- 

To further explore the VL jamming hypothesis, we 
performed a series of Small Angle Neutron Scattering 
(SANS) experiments. By preparing a metastable VL and 
then inducing vortex motion by small changes in the mag- 
netic field, we resolved coexisting metastable and ground 
state phases and obtained a quantitative measurement 
of the VL transition. Our results provide conclusive evi- 
dence that metastability in the VL phases of MgB2 is not 
a result of pinning to defects and lend further credibility 



to a jamming scenario. 

The SANS experiments were performed on the Dll 
beam line at the Institut Laue-Langevin (ILL) and on the 
CG2 General Purpose SANS beam line at the High Flux 
Isotope Reactor (HFIR) at Oak Ridge National Labo- 
ratory. To achieve diffraction, the sample and magnet 
were rotated and/or tilted together in order to satisfy 
the Bragg condition for the VL planes. To resolve closely 
located VL Bragg reflections, very tight collimations of 
the neutron beam were used (D22: 0.03°, HFIR: 0.07° 
FWHM). Measurements were performed using the same 
200/ig MgB2 single crystal as in previous SANS exper- 
iments l8[ . The sample was grown using isotopically 
enriched 11 B to decrease neutron absorption and had 
a critical temperature T c = 38 K and upper critical field 
H C 2 = 3.1 T. Measurements were performed at 2 K and 
0.5 T applied parallel to the c axis. 

The ground state (GS) VL phase diagram for MgB2, 
shown in Fig. QJa), consists of three different hexag- 
onal configurations. In the F phase, VL Bragg peaks 
are aligned along the crystal a axis, and in the I phase 
along the crystal a* axis. In the intermediate L phase, 
the VL rotates continuously from the a to a* orientation 
[H, [Toj| . Examples of diffraction patterns for the F and L 
VL phases are shown in Figs. 0Jb) and (c) respectively. 
The degeneracy in rotation direction of the VL domains 
in the L phase results in 12 diffraction peaks, unlike the 
6 peaks observed in the F and I phases. For a given field 
and temperature, the GS VL was reached by applying 
a damped oscillation of the magnetic field with a initial 
amplitude ~50 mT. A metastable (MS) VL is prepared 
by heating or cooling across the F-L or L-I phase bound- 
aries. The MS VL diffraction pattern in Fig. 0Jb) was 
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FIG. 1. MgB2 vortex lattice phases for H \\ c. (a) Ground 
state phase diagram. Diffraction patterns in (b) and (c) were 
measured at 2 K and 0.5 T, indicated by the red star. The 
metastable VL (b) is formed by cooling across the F-L phase 
boundary as shown by the arrow. A 50 mT damped field 
oscillation drives the VL to the ground state configuration 
(c). Crystalline axes are shown in (c). 



obtained by preparing a GS VL at 17 K and 0.5 T, and 
then cooling to 2 K as indicated by the arrow in Fig.QJa). 
A field oscillation at 2 K yielded the GS VL diffraction 
pattern in Fig. DJc). 

By subjecting the MS VL to small changes in the 
magnetic field, a gradual transition to the GS is ob- 
served. Figure (2fa) shows the as prepared MS VL at 
2 K and 0.5 T. After reducing the applied magnetic field 
by 22 mT, VL Bragg peaks corresponding to both the 
MS and GS phases are present, as seen in Fig.[2fb). The 
tight collimation used in the experiment allows a clear 
resolution of the individual peaks and shows that MS and 
GS VL domains coexist within the sample. In Fig. [2fc) 
the magnetic field has been further reduced by a total 
of 65 mT, driving the VL to the GS within the entire 
sample. 

Measurements of the relative intensity of the diffrac- 
tion peaks makes it possible to determine the fraction of 
the VL in the GS and MS states. Figures [^d-f) show the 
azimuthal intensity distribution corresponding to panels 
(a-c). The peaks were fitted using three gaussians with 
identical widths and peak centers fixed at 259.9°, 265.8°, 
and 271.8° respectively. The fitted areas under the three 
peaks (Aqsi, Ams, ^GS2) provide a quantitative measure 
of the population of MS and GS VL phases in the sample 
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FIG. 2. Resolving the MS to GS transition, (a-c) Diffraction 
patterns at 2 K: As prepared MS VL state at 0.5 T (a); Co- 
existence of MS and GS obtained after decreasing the field by 
22 mT (b); GS VL obtained after decreasing field by 65 mT 
(c). (d-f) Azimuthal intensity distribution corresponding to 
the diffraction patterns. Solid lines are fits to the data as 
described in the text. 



and allow the calculation the relative volume fractions: 

' ^MS + ^GS2) 



/ms = A M s I (^GSl 
/gs = 1 - /ms- 



(1) 

(2) 



Figure [31(a) shows the evolution of the MS and GS VL 
volume fraction as the vortex lattice is driven from the 
metastable to the ground state. Initially, an MS VL was 
prepared at 0.5 T and 2 K, followed by decreases in the 
applied field in steps of 2 to 5 mT. The VL was remea- 
sured after every change of the magnetic field, and /ms 
and /gs were calculated using Eqns. (pQ) and (|2]). The 
metastable VL volume fraction smoothly decreased as 
the applied field was decreased. 

An straightforward explanation for VL metastability is 
that it may be due to vortex pinning. In this picture the 
VL transition to the ground state is due to vortex mo- 
tion induced by the decreasing field, as the front of vortex 
motion propagates from the edges of the crystal toward 
the center. The red line in Fig. [31(a) is calculated using a 
simple two dimensional Bean pinning model with a char- 
acteristic field hqH* = 50 mT Although this pro- 
vides an overall good fit, there are noticeable deviations 
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FIG. 3. VL transition from MS to GS. (a) Relative inten- 
sity of the MS and GS VL Bragg peaks measured at 2 K as 
the applied field was decreased from 0.5 T. The black line 
is an exponential fit to the data with a characteristic field 
H* = 18 mT, and the red line is a Bean model fit with a 
critical field H* = 50 mT. (b) Solid symbols correspond to an 
initial decrease of 14 mT followed by an increase of the ap- 
plied magnetic field. Open symbols same as in (a), (c) Same 
as (b) but with solid symbols reflected about the reversal field 
of 0.486 T. 



in the field ranges 0.43-0.46 T and 0.47-0.49 T. In con- 
trast, a simple exponential function with a characteristic 
field of 18 mT (black line) provides a better fit through- 
out the entire field range. In addition, the characteristic 
field predicted by the Bean model n^H* = {1qJ c D/2, 



where J c ~ 10 6 A/m [12|, ll3| and the sample diameter 
D ~ 1 mm, yields a fi^H* ~ 1 m T, a value more than an 
order of magnitude smaller than that suggested by the 
data in Fig. [3fa). 



To further investigate whether the VL metastability 
could be due to vortex pinning, a second sequence of mea- 
surements was performed using a field reversal. Within 
the Bean picture, a field reversal will induce a second, 
inward moving, flux flow front, and the further transi- 
tion to the VL ground state should not occur until this 
front reaches the metastable portion in the center of the 
sample. A MS VL was prepared and the applied field de- 
creased by 14 mT which rendered the VL in a state with 
/ms = 53%. The field was then increased in 2 mT incre- 
ments, causing vortices to reenter the sample from the 
edges. The results from the field reversal measurements 
are shown in Fig. [3jb). Contrary to the Bean model pre- 
diction, no plateau was observed in the population of the 
MS VL. Rather, when the data is reflected about the re- 
versal field of 0.486 T, it coincides within error bars with 
the results from the decreasing magnetic field measure- 
ments, as shown in Fig.[3fc). This shows that any change 
in the applied field on the order of a few millitesla (con- 
sistent with our estimate of H*) is sufficient to perturb 
some fraction of the metastable VL domains and cause a 
transition to the ground state. 

It is possible that the vortex motion is not well de- 
scribed by the simple Bean model. Magneto-optical mea- 
surements on thin films of various superconductors in- 
cluding MgB2 have shown that, in some cases, vortices 
enter through dendritic avalanches which extend towards 
the center of the sample and then gradually fill the entire 
volume 



In a similar scenario, one could imagine a 
situation with vortex motion confined to dendrites within 
which the VL has reoriented to the ground state. How- 
ever, it is important to note that the observed dendritic 
instability is associated with the initial flux entry into 
the thin film, whereas our measurements are carried out 
with a uniform vortex density throughout the sample. 

To definitively resolve whether the VL metastability is 
due to pinning, we considered the vortex density corre- 
sponding to the measurements presented in Fig. [3j For 
a hexagonal vortex lattice, the VL scattering vector q 
depends on the magnetic induction B as 
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The scattering vector, and thus B, can be determined 
experimentally from the peak positions in the VL diffrac- 
tion patterns. Figure 2] shows the magnetic induction as 
a function of applied field corresponding to the data in 
Fig. [3^b). Within the scatter in the data, the magnetic 
induction for both the MS and GS VL domains is found to 
be proportional to the applied field. The small deviation 
(1.4%) of B I hqH from unity is within the experimental 
uncertainty. The error on the determination of B ex- 
ceeds our estimate of fioH* ~ 1 mT, and a radial peak 
broadening due to the field variation within the sample 
is therefore not observed. In contrast, a VL pinned by 
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FIG. 4. Magnetic induction calculated from the measured 
VL scattering vector as a function of applied magnetic field, 
corresponding to the data in Fig.[3](b). The dashed line shows 
a linear fit to the data and the solid line is B — jiqH. 




main jamming represents a novel type of collective vor- 
tex behavior, distinct from the jamming of individual 
vortices observed in artificially engineered pinning po- 
tentials in vortex ratchets and similar devices [3-6] . Do- 
main jamming would most likely also be reflected in the 
dynamical properties of the VL. Analogous to the jam- 
ming observed in granular materials one might expect 
the emergence of power law behavior 16l-[l8|. and it is 
possible that the slow relaxation of the VL found in the 
non-centrosymmetric superconductor Li2PtaB is a con- 
sequence of VL domain jamming [l9| . 

In summary, we have made the first detailed study of 
the metastable VL in MgE>2 while it was gradually driven 
to the ground state by small decreases in the applied 
magnetic field. Our measurements provide definitive ev- 
idence that the metastability can not be ascribed to vor- 
tex pinning. Instead, we propose that the metastability 
in MgE>2 results from the jamming of counter rotating 
VL domains. Further work is required to explore this 
hypothesis. 
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actor was sponsored by the Scientific User Facilities Di- 
vision, Office of Basic Energy Sciences, U. S. Department 
of Energy. 



FIG. 5. Schematic illustrating the proposed jamming of 
counter-rotated VL domains. For the center domain, rotation 
to the ground state orientation is prevented, as it would lead 
to overlapping domains in the shaded regions. 



defects would have a fixed vortex density, independent of 
changes in the applied field. This unequivocally proves 
that a mechanism other than pinning is responsible for 
the VL metastability in MgB2. 

We propose that the metastability is due to a jamming 
of counter rotated VL domains which prevents a rota- 
tion to the ground state orientation, as shown schemat- 
ically in Fig. [5] In this scenario the VL domain bound- 
aries, arising from the two degenerate orientations of the 
ground state (L phase), must be robust enough to sup- 
port a jammed state and ensure the longevity of the 
metastable VL phases. Furthermore, the domain bound- 
aries must persist even as the vortex lattice expands or 
is compressed as a result of the changing magnetic field. 
Theoretical studies of VL domain boundaries have, to 
date, been limited [Hj, and further work is required to 
determine the feasibility of the proposed VL domain jam- 
ming. It is important to note that the proposed VL do- 
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